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Abstract 

These are the proceedings of the workshop "HadAtom99" , held at the Institut fiir Theoretische Physik, 
Universitat Bern, October 14 - 15, 1999. The main topics discussed at the workshop were the physics of 
hadronic atoms and in this context recent results from Chiral Perturbation Theory. Included here are the 
program, the list of participants and a short contribution from the speakers. 
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1 Introduction 



There has been a growth of interest in hadronic atoms in recent years. On the experimental side, energy levels 
of hadronic atoms as well as their decay into neutral hadrons have been or will be measured with high precision: 
Pionic hydrogen and n^d atoms at PSI, pionium at CERN (DIRAC), and K^p as well as K~d atoms at 
DAFNE (DEAR). On the theoretical side, effective lagrangian techniques have shown to be very efficient also 
here: Bound state calculations - in particular the determination of energy levels and decay widths - can now 
be performed with surprising ease. 

For these reasons, Leonid Nemenov suggested to organize a workshop on this subject at the Institut fiir 
Theoretische Physik, Universitat Bern, similarly to the one organized at Dubna (Russia) in May 1998 on the 
same subject. The talks were devoted to recent experimental and theoretical progress in the investigations of 
hadronic atoms. 

The meeting had 48 participants whose names, institutes and e-mail addresses are listed below. Results were 
presented in 25 minutes talks by 26 physicists. A one page abstract of these presentations as well as a list of 
the most relevant references has been provided by most of the speakers. These contributions are collected here. 
Below we display the list of participants with their e-mail, then follows the program of the workshop. Finally, 
we provide the individual excerpts submitted by the speakers 0. 

We would like to thank all participants for their effort to travel to Bern and for making this an exciting and 
lively meeting, and our secretaries, Ruth Bestgen and Ottilia Hanni, for their excellent support. Last but not 
least, we thank our colleagues from the organizing committee (Heinrich Leutwyler, Valery Lyubovitskij , Leonid 
Nemenov, Hagop Sazdjian, and Dirk Trautmann) for their invaluable help in structuring the meeting. 



Bern, November 1999 



Jiirg Gasser, Akaki Rusetsky and Jiirg Schacher 



^ This kind of miniproceedings had been collected and posted on the archive in recent years by Hans Bijnens and Ulf Meifiner 
at several occasions, see e.g. J. Bijnens and Ulf-G. Meifirinr, Chiral Effective Theories, 205. WE-Heraeus Seminar, Bad Honnef, 



November 30 - December 4, 1998, Miniproceedings, hcp-ph/990138] . We found this an excellent way to summarize the event: 
does not demand excessive labor from the authors, yet it provides a lot of useful information. Here we follow this idea. 



It 
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3 HadAtom99: Program 



Thursday, October 14 



10:00 - 10:05 J. Gasser (Bern) 

10:05 - 10:30 J. Schacher (Bern) 

10:30 - 10:55 R. Lednicky (Prague) 

10:55 - 11:20 A.V. Tarasov (Dubna) 

11:20 - 11:45 Coffee 

11:45 - 12:10 F. Ravndal (Oslo) 

12:10 - 12:35 J. Soto (Barcelona) 

12:35 - 13:00 H. Sazdjian (Orsay) 

13:00 - 14:30 Lunch 

14:30 - 14:55 V.E. Lyubovltskij (Dubna) 

14:55 - 15:20 A. Rusetsky (Bern) 

15:20 - 15:45 L.G. Afanasyev (Dubna) 

15:45 - 16:15 Coffee 

D. Trautmann (Basel) 

16:15 - 17:15 Th. Heim (Basel) 
K. Mencken (Basel) 

17:15 - 17:35 Coffee 

17:35 - 18:00 G. Baur (FZ Julich) 

18:00 - 18:25 J. Cugnon (Liege) 



Welcome 

The DIRAC experiment at CERN 

Finite size effects on tt+tt" production in 
continuous and discrete spectrum 

Analytical approach to calculation 
of the TT+TT" atom production 
rate in different quantum states 



Corrections to the pionium lifetime 

Effective field theory approach to pionium 
Pionium lifetime in Generalized Chiral 
Perturbation Theory 



TT+TT atom in QCD 

Spectrum and decays of hadronic atoms 
Calculation of the breakup probability of 
7r+7r~ atom in a target with a high accuracy 



Group report 
about pionium 
interaction with matter 



Another exotic relativistic atom: antihydrogen 

Multiphoton exchange effects in 
pionium-matter interaction 
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Friday, October 15 



9:00 - 


9:25 


L.L. Nemenov (CERN-Dubna) 


Development in the theoretical dimesoatom 
description required by DIRAC experiment 


9:25 - 


9:50 


C. Guaraldo (LNF-INFN) 


The DEAR experiment at DA$NE 


9:50 - 


10:15 


L. Simons (PSI) 


Pionic hydrogen: status and outlook 


10:15 - 


■ 10:40 


Coffee 




10:40 - 


■ 11:05 


U.-G. MeiBner (FZ Jiilich) 


Pion-kaon scattering 


11:05 - 


■ 11:30 


H. Leutwyler (Bern) 


Discussion: pion-kaon scattering 


11:30 - 


11:45 


Coffee 




11:45 - 


■ 12:10 


P. de Simone (LNF-INFN) 


Kii decays at DA$NE 


12:10 - 


12:35 


S. Pislak (BNL) 


A new measurement of the 

— > 7r+7r~e+i/ decay 


12:35 - 


■ 13:00 


G. Amoros (Lund) 


Ku decays: a theoretical discussion 


13:00 - 


■ 14:30 


Lunch 




14:30 - 


■ 14:55 


J. Stern (Orsay) 


Chiral phase transitions at zero temperature 


14:55 - 


■ 15:20 


L. Girlanda (Orsay) 


The two-flavor chiral condensate from 
low-energy tttt phase-shifts 


15:20 - 


15:45 


M. Knecht (Marseille) 


Radiative corrections in Ki4 decays 


15:45 - 


■ 16:15 


Coffee 




16:15 - 


■ 16:40 


G. Colangelo (Zurich) 


Numerical solutions of Roy equations 


16:40 




J. Gasser (Bern) 


Closing remarks 
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The DIRAC Experiment at CERN 



DIRAC Collaboration 
J. Schacher 

Laboratorium fur Hochenergiephysik, Universitdt Bern, CH-3012 Bern, Switzerland 



The experiment DIRAC [1], a magnetic double arm spectrometer, aims to measure the tt+tt" atom lifetime 
with 10% precision, using the high intensity 24 GeV/c proton beam of the CERN Proton Synchrotron. Since 
the value of this lifetime of order 10~^^ s is dictated by strong interaction at low energy, a precise measurement 
of this quantity enables to study characteristic pion parameters in detail and to submit predictions of QCD to 
a severe check. 

After tmiing the' primary proton beam as well as all detector components in detail, the experiment started to 
accumulate data this summer. Preliminary results concerning quality and reliability of the setup are presented 
in the following two figures: 
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Figure 1: Correlated pairs of 7r+ and tt with 
momenta piab smaller than 4.5 GeV/c and 
Qt < 4 MeV/c. 
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Figure 2: Invariant pair mass of tt and p with 
momentum piab larger than 3 GeV/c. 



In Fig. 1, the distribution of tt+tt pairs in dependence oi Ql shows a noticeable peak in the yield of tt+tt 
pairs with momenta Q less than 5 McV/c, due to Coulomb attraction in the final state. The longitudinal 
component Ql is defined as the projection of the relative momentum Q in the tt+tt^ system on the direction of 
the TT+TT" lab momentum. The mean value of the peak, (Ql) = 0.095 MeV/c, is small and hence clear evidence 
for good setup alignment. 

As emphasized in Ref. 1, the resolution in Ql and Qt plays a key role in the data analysis. In order to keep 
track of the relative momentum resolution, the invariant mass distribution of p and 7r~ — products of the A 
decay was carefully investigated. The position of the mass peak in Fig. 2 corresponds to A particles, recorded 
and reconstructed in the magnetic spectrometer. The width of the A peak, mainly given by the momentum 
resolution, is found to be less than 1 MeV/(? (see Fig. 2). Thus, the DIRAC collaboration is confident, that 
the magnetic field and the location of the track detectors are known with adequate precision. 



[1] B. Adeva et al., Proposal to the SPSLC, CERN/SPSLC 95-1 (1995). 
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Finite— size effects on tt+tt production in continuous and discrete 

spectrum 

R. Lednicky 

Institute of Physics, Na Slovance 2, 18221 Prague 8, Czech Repubhc 

The determination, on a per cent level accuracy, of the breakup probability P^r ~ /Na of the tt+tt" 
atoms produced by a high energy beam in a target is of principle importance for a precise lifetime measurement 
of these atoms in the experiment DIRAC [1] currently running at CERN. Clearly, the breakup probability is a 
unique function of the target geometry and material, the Lorentz factor and the ground state lifetime. 

While the number of the breakup atoms N^" is measured directly as an excess of the tt^tt^ pairs with very 
low relative momentum, the number iV^ of the produced tt+tt" atoms is calculated based on the Migdal- Watson 
theory of final state interaction in discrete [2,3] and continuous [4] spectrum. This calculation is sensitive to 
the space-time extent of the pion production region, characterized by the distance r* between the 7r+ and tt~ 
production points in their c.m.s.. In ref. [2], the r* -dependence was treated in an approximate way, dividing 
the pion sources into short-lived and long-lived ones. It was assumed that r* = for pion pairs arising solely 
from the short-lived sources and characterized by the distances r* much smaller than the Bohr radius of the 
7r+7r~ system: \a\ = 387.5 fm, otherwise r* = oo. 

The finite size correction to such calculated number of free tt+tt" pairs arises mainly from the region of 
small relative momenta compared with 1/r*. It is determined by the three dimensionless combinations r* /a, 
fo/r* and fo/a of r*, a and the s-wave 7r+7r^-scattering length: /o = ■|(2no + CL2) ~ 0.23 fm. Typically (r*) 
10 fm so that the correction is dominated by the strong interaction effect and can amount up to ~ 10%. 

Due to a small binding energy e ~ (2/xa^)~^, the finite-size correction to the discrete spectrum at r* <C \a\ 
is nearly the same as that to the continuous spectrum at zero energy. Since Na is actually determined by a 
ratio of the pions produced in discrete and continuous spectrum, this correction, up to a fraction of r* /a, would 
cancel out in the breakup probability provided we could measure the number of free tt+tt" pairs in the region of 
very small relative momenta, much less than 1/r* [5,6]. With the increasing relative momentum the correction 
to the production of the free pairs decreases, the decrease being faster for pions emitted at larger distances r* . 
Particularly important is the case when one of the two pions comes from the w-meson decay and r* ^ 30 fm. 
As a result, in DIRAC conditions, we expect the net finite size correction in the breakup probability of the 
order of per cent. It leads to several per cent effect in the measured lifetime and, in first approximation, can be 
neglected compared with the expected statistical error of ~ 10%. 

We discuss how to diminish the systematic error related to the finite size effects, using the correlation data 
on identical charged pions (containing the information about the distances r* between the pion production 
points in the same experiment) together with the complete phase-space simulations within transport models. 
We also show that the usual equal-time approximation, the neglect of the space-time coherence as well as of 
the transition tt+tt" 7r°7r° and of the residual nucleus charge introduce negligible systematic errors of the 
order of per mill in the measured lifetime. 

[1] B. Adeva et al.: Lifetime measurement of tt+tt" atoms to test low energy QCD predictions. 
Proposal to the SPSLC, CERN/SPSLC 95-1, updated 10 November 1995. 
[2] L. Nemenov: Yad. Fiz. 41 (1985) 980. 

[3] V.L. Lyuboshitz: Yad. Fiz. 48 (1988) 1501 {Sov. J. Mud. Phys. 48 (1988) 956). 

[4] R. Lednicky, V.L. Lyuboshitz: Yad. Fiz. 35 (1982) 1316 (Sov. J. Nucl. Phys. 35 (1982) 770). 

[5] R. Lednicky: On the breakup probability of tt+tt" atoms, DIRAC meeting, February, 1998; 

Finite-size effects on tt+tt" production in continuous and discrete spectrum, DIRAC meeting, February, 1999. 

[6] L. Afanasyev, O. Voskresenskaya: Phys. Lett. B 453 (1999) 302. 
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Analytical approach to calculation of the tt+tt atom production rate 

in different quantum states 

A.Tarasov ^'^ and O.Voskresenskaya^'^ 
-^Institut fiir Theoretischc Physik, Universitat Regensburg, D 93040 Regensburg, Germany 
^Max-Plank Institut fiir Kernphysik, Postfach 103980, D-69029 Heidelberg, Germany 
^Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 

A method of measurement of the 7r+7r~-atom lifetime proposed by L.L.Nemenov [1] and which is under 

realization in the experiment DIRAC [2] at CERN Proton Synchrotron is essentially based on the assumption 
that the ratio R{nS/pS) of the 7r"'"7r~-atom production rate in nS'-states in hadron-nuclear interactions at high 
energy to the production rate of so the called "Coulomb" tt+tt" pairs (i.e. free tt+tt" pairs produced at a 
distance less or of the order of the 7r"'"7r~-atom Bohr radius and having the relative momentum of an order of 
the Bohr momentum) can be calculated with accuracy better than 1%. Here we consider the accuracy of this 
ratio. 

According to [1,3], the ratio can be written as: 



where ipnsi'r) is the 7r"'"7r~-atom wave function with the principal quantum number n and zero orbital one (ns- 
states); ipps{r) is the wave function of tt+tt" pair with the relative momentum p and with zero orbital momentum 
(s-wave); M(f) is the amplitude of 7r"'"7r~-systcm production at the relative distance f. 

In ref.[4] it is shown that in the first order of the perturbation theory over the strong interaction R{nS/pS) 
is expressed via the a squared ratio of the well-known Coulomb wave functions. 

Numerical calculations in the resent paper [5] show that the same ratio between discrete states is also 
expressed via the a squared ratio of the Coulomb wave functions at origin. 

In this paper we analyse the influence of the strong interaction on the ratio R(nS/pS) analytically using 
a square well approximation for the strong interaction potential. It is shown that the result of ref.[4] is valid 
with accuracy of an order of ratio of the tt"*"?:" -scattering length to the 7r"'"7r~-atom Bohr radius which is 
approximately equal to 10~^. Thus, both perturbative and nonperturbative calculations of R{nS/pS) confirm 
that Nemenov's result [1] for this quantity is valid with the required accuracy. 

[1] L.L. Nemenov, Yad.Fiz. 41 (1985) 980. 

[2] B. Adeva et al., Lifetime measurement o/tt+tt" atom,s to test low energy QCD predictions, (Proposal to the 
SPSLC, CERN/SPSLC 95-1, SPSLC/P 284, Geneva 1995). 

[3] L.G. Afanasyev, 0.0. Voskresenskaya and V.V. Yazkov, JINR Communication, Pl-97-306, Dunba 1997. 
[4] L.G. Afanasyev and 0.0. Voskresenskaya, Physics Letters B 453 (1999) 302. [5] I. Amirkhanov, I. Puzynin, 
A. Tarasov, O. Voskresenskaya and O. Zeinalova Physics Letters B 452 (1999) 1. 
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Corrections to the pionium lifetime 



Finn Ravndal 

Institute of Physics, University of Oslo, N-0316 Oslo, Norway. 



The charged pions with mass m+ in the IS ground state of pionium are bound in the Coulomb wavefunction 
\E'(r) with characteristic momentum 7 = Its dominant decay mode is into two neutral pions with mass 

Too- Since the mass difference Am ~ m+ — mg is small, both the pions in the initial state of the hadronic 
atom and in the final state of the decay are non-relativistic. This was first emphasized by Labelle and Buckley 
[1] who showed that the leading term of the decay rate F can be obtained by calculating the energy shift Ai? 
within non-relativistic effective field theory using the relation F = — 2ImAi?. The corresponding Lagrangian 
was then extended by Kong and Ravndal [2] to include higher order interactions which allow corrections to the 
dominant decay rate to be systematically calculated. 

In order to determine the unknown coupling constants of the effective Lagrangian, one must match the non- 
relativistic scattering amplitude for 7r+ + tt^ ^ tt'^ + 7r° from this theory to that which follows from relativistic, 
chiral perturbation theory at the same energy. At low momenta this latter one is defined to be 

TflXp) = 327r(a + 6pV"^l) (1) 

in the center-of-mass frame where the charged pions have the momentum p. The scattering length a and the 
slope parameter h include both higher order chiral corrections and isospin- violating effects from the quark mass 
difference vriu — md and short-range electromagnetic effects. In the matching one must also include the effects 
coming from the difference in normalization of relativistic and non-relativistic states as done correctly by Gall, 
Gasser, Lyubovitskij and Rusetsky [3]. 

At threshold the momentum p — and the full scattering amplitude with long-range Coulomb interactions 
removed is then given by just this scattering length. It can be written in terms of the more conventional isospin- 
symmetric scattering lengths ao and 02 in the isospin / = and isospin 1 = 2 channels as a = qq — a2 + Aa 
where Aa includes these symmetry-breaking effects. 

Calculating now the transition rate and multiplying by the probability |\I'(0)p = T'^/tt that the pions are 
found at the same point, one obtains the leading order result for the decay rate 

Fo = la^'^V^Ammoa^ (2) 
9 m+ 

Relativistic effects in the final state can be obtained from a corresponding term of the effective Lagrangian [4] 
and follows also directly from the phase space factor 

5{2m+ - 2Eo) (3) 

where Eo = rriQ + A;^/2mo — k^/SrriQ is the energy of the neutral pions when we ignore the small binding energy. 
To this order in Am/m the decay rate is thus found to be 



A 



F = ^aV2Ammo i^l + (4) 

This is consistent with the more general result of Gall et al. [3] and has also been found by Eiras and Soto [5] . 

Remaining corrections can now be obtained by including long-range Coulomb interactions and rescattering 
effects. These have been worked out by Gall et al. [3] and will also follow from the effective theory of Eiras and 
Soto [5]. The above non-relativistic approach can be similarly applied to other hadronic atoms. 



[1] P. Labelle and K. Buckley, ^iep-ph/9804201 



[2] X. Kong and F. Ravndal, Phys. Rev. D59, 014031 (1999). 
[3] A. Gall, J. Gasser, V.E. L yubovitskij and A. Rusetsky, Phys. Lett. B462, 335 (1999). 



[4] X. Kong and F. Rav ndal, |hep-ph/990553£ 



[5] D. Eiras and J. Soto, |hep-ph/990554£ 
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Effective field theory approach to pionium 

Joan Soto 

Departament d'Estructura i Constituents de la Materia, Universitat de Barcelona, Diagonal 647, 08028 

Barcelona, Catalonia, Spain 



I shall briefly report on the model-independent approach to pionium which has been presented in [1] . 

It basically consist of identifying the relevant energy and momentum scales of the tt+tt" atom and sequentially 
integrate them out until we reach the lowest relevant scale, namely the binding energy. This is carried out by 
introducing a series of non-relativistic effective field theories and by requiring them to be equivalent at the 
desired order of accuracy. Since pionium is an electromagnetic bound state, it contains at least three dynamical 
scales: the mass of the charged pions m ~ lAOMeV, the typical relative momentum in the bound state 
Tna/2 ^ 0.5MeV, and the binding energy ^ 2keV [2]. In addition pionium decays mainly to two neutral 

pions through the strong interactions, which brings in an additional energy scale Am ^ 5MeV, the difference 
between the charge and neutral pion masses, and its associated three momentum s = \/2mAm ~ AOMeV. 
In view of these numerical values two important observations are in order: (i) all the scales are considerably 
smaller than m ~ 140 MeV, which implies that a non-relativistic approach should be appropriated, and (ii) all 
the scales, including m, are considerably smaller than the typical hadronic scale, say the rho mass, and hence 
all the necessary information is contained in the Chiral Lagrangian (xL) coupled to electromagnetism. Since 
at these energies not only electromagnetism but also the strong interactions are amenable to a perturbative 
treatment (in a and in Chiral Perturbation Theory (xPT) respectively), pionium can be addressed in a model 
independent way. 

Let us then start with the SU{2) x SU{2) xL coupled to electromagnetism. After integrating out the 
scale m for pion pairs near threshold one should obtain what may be called a Non-Relativistic xL (NRxL). 
The degrees of freedom of NRxL are photons and non-relativistic charged and neutral pion fields. NRxL is 
analogous to Non-Relativistic QED [2] , but contains non-relativistic pseudoscalar fields instead of Pauli spinor 
fields and includes the effect of the strong interactions. We present the general form of the NRxL including 
isospin breaking terms. Since Am 3> ma/2 , mQ;^/4 we can also integrate out this scale and its associated three 
momentum s (neutral pions). We obtain a second Non-Rclativistic xL which may be abbreviated as NRxL'. 
Its degrees of freedom are photons and non-relativistic charged pions. The introduction of NRxL' avoids having 
to solve eventually a coupled channel problem. Finally, integrating out photons of energy or momentum of the 
order ma/2 we obtain what may be called potential NRxL' (pNRxL'), in analogy of potential NRQED [3]. At 
the order we are interested in, pNRxL' contains non-relativistic charged pions interacting through a potential 
(electromagnetic and strong), and it is totally equivalent to a suitable quantum mechanical Hamiltonian. 

A second order quantum mechanical perturbation theory calculation in pNRxL' allows us to obtain the 
decay width and the spectrum at next to leading order in Am/m, ma'^/4Am and a. In order to match this 
precision a yet-to-be-done two loop calculation in the xL including photons is necessary. 

The techniques described above may also be useful for other hadronic atoms. 

Related work using non-relativistic effective field theories can be found in [4]. 



[1] D. Eiras and J. Soto, |hep-ph/9905"54^ 



[2] W.E. CasweU and CP. Lepage, Phys. Lett. B167 (1986) 437. 



X. Kong and F Ravndal, |hep-ph/9805357 



B. R. Holstein, aucl-th/9901041 



[3] A. Pineda and J. Soto, Ph ys. Rev. D59 (199 9) 016005, |hep-ph/9805424 
[4] P. Labelle and K. Buckle y ^iep-ph/980420"ll 



A. GaU, J. Gasser, V. E. Lyubovitskij and A. Rusetsky, Phys. Lett. B462 (1999) 335. 
X. Kong and F. Ravndal, |hep-ph/9905539 . 
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Pionium lifetime in generalized chiral perturbation theory 

H. Sazdjian 

Groupe de Physique Theorique, Institut de Physique Nucleaire, 
Universite Paris XI, 91406 Orsay Cedex, Prance 

The relationship between the pionium hfetime and the tttt scattering lengths, including the sizable electro- 
magnetic corrections, is analyzed in the framework of generalized chiral perturbation theory, in which the quark 
condensate value is left as a free parameter. The variation curve of the lifetime as a function of the combination 
(oq — Oq) of the S-w&ve scattering lengths is presented. 
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TT+TT Atom in QCD 

J. Gasser\ V.E. Lyubovitskij ^-^, and A.G. Rusetsky^-^'^ 
^ Institute for Theoretical Physics, University of Bern, Sidlerstrasse 5, CH-3012, Bern, Switzerland 
^ Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia 
^ Department of Physics, Tomsk State University, 634050 Tomsk, Russia 
^ HEPI, TbiUsi State University, 380086 TbiUsi, Georgia 



The DIRAC experiment at CERN [1] aims to measure the lifetime of the tt+tt" atom in its ground state 

with high precision. The atom predominantly decays into two neutral pious. The measurement of the quantity 
r2^o will allow one to determine the difference |ao — 02] of the strong S'-wave tttt scattering lengths with isospin 
/ = 0, 2 and to check the predictions for this quantity obtained in the standard version of ChPT [2] , and to 
investigate the nature of spontaneous chiral symmetry breaking in QCD [3]. In order to perform this programme, 
the theoretical expression for the width must of course be known with a precision that matches the accuracy of 
the lifetime measurement of DIRAC. 

We present a general expression for the corresponding decay width in the framework of QCD (including 
photons). This expression is proportional to the square of the relativistic 7r"'"7r~ 7r°7r° scatterting amplitude 
A, contains all terms at leading and next-to-leading order in isospin breaking, and is valid to all orders in the 
chiral expansion [4] : 

r2.o = ^aV^'(l + ^), (1) 

where p* = (M^+ — M^o — c? M'^+ / 4)^/'^ . The quantities A and K are expanded in powers of a and — m„. 
We count a and (m^ — m„)^ as small parameters of order 5: 

A = ao - 02 + /ii (m^ - m„)^ + /i2 a + o(5), 

K = h{md-muf + f2a\na + ha + o{5), (2) 

where the scattering lengths ao, 0,2 and the coefficients hi and fi are evaluated in the isospin symmetry limit, 
i.e. at m„ = = m and = = 139.57 MeV. The quantities hi and parameterize the corrections 
to the leading-order Deser's et al. type [5] formula. For the factor K we get the exact expression (without any 
chiral expansion) 

K = (ao + 2a2f + ^ (1 - Ina) (2ao + 02) + o{5). (3) 

7r+ 

Therefore, the formula (1) relegates the problem of the calculation of the tt+tt" atom lifetime to the evaluation 
of the physical on-mass-shell scattering amplitude for the process tt+tt" — > 7r°7r° to any desired order in the 
chiral expansion. 

Next we apply Chiral Perturbation Theory at one loop to analyse the formula for the tt+tt" atom lifetime 
(1). We obtain the analytic expressions for the next-to-leading order 0{a) and 0{md — niuf corrections and 

pin them down numerically. 

[1] B. Adcva et al, CERN proposal CERN/SPSLC 95-1 (1995). 

[2] J. Gasscr and H. Leutwyler, Phys. Lett. 125B (1983) 325; J. Bijncns, G. Colangelo, G. Ecker, J. Gasser, 
and M. E. Sainio, Phys. Lett. B374 (1996) 210. 

[3] M. Knccht, B. Moussallam, J. Stern, and N. H. Fuchs, Nucl. Phys. B 457 (1995) 513; ibid. B 471 (1996) 
445. See also Girlanda's talk at this meeting. 

[4] A. Gall, J. Gasser, V. E. Lyubovitskij, and A. Rusetsky, Phys. Lett. B 462 (1999) 335. 
[5] S. Dcscr, M. L. Goldberger, K. Baumann, and W. Thirring. Phys. Rev. 96 (1954) 774; J. L. Uretsky and 
T. R. Palfrey, Jr., Phys. Rev. 121 (1961) 1798; S. M. Bilenky, Van Kheu Nguyen, L. L. Nemenov, and 
F. G. Tkebuchava, Yad. Fiz. 10 (1969) 812 . 
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Spectrum and decays of hadronic atoms 



J. Gasser\ V.E. Lyubovitskij^'^, and A.G. Rusetsky ^-^-''' 
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^ Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia 
^ Department of Physics, Tomsk State University, 634050 Tomsk, Russia 
* HEPI, TbiUsi State University, 380086 TbiUsi, Georgia 

A general survey of the theoretical approaches to the description of hadronic atom characteristics (spectrum 
and decay width) is given. We confront nonrelativistic effective Lagrangian approach, as well as other field- 
theoretical approaches based on relativistic equations, with the potential scattering theory formalism for the 
hadronic atom problem. This comparison is essential in the view of several ongoing or planned experiments 
(DIRAC at CERN. PSL DEAR at LNF-INFN, KEK, Uppsala) that are aimed at the high-precision measurement 
of hadronic atom characteristics. These measurements will enable one to directly extract the strong scattering 
amplitudes from the data, and thus contribute a valuable piece of information to our knowledge of the dynamics 
of strong interactions at low energy. 

The analysis of the data in the above-mentioned experiments are carried out with the use of the relations 
given in Ref. [1]. Up to the corrections in the next order in isospin breaking, these relate the energy- level 
shift of the hadronic atom llE and its partial decay width into the channel with neutral isotopic partners, 
to the certain isotopic combinations of the strong scattering lengths. However, in order to fully exploit the 
high-precision data available from recent and future experiments on hadronic atoms, one has to consistently 
evaluate the corrections to the relations given in Ref. [1]. Using the powerful technique based on the effective 
nonrelativistic Lagrangians, it is demonstrated that in the next-to-leading order in isospin breaking, and in all 
orders in chiral expansion, the partial decay width of tt+tt" atom into is given by [2] 

r^^o = la^Ml,-Ml.-\a^Ml,f/^A'{l + K) 

K = g^(ao + 2a2)2 + ^(l-lna)(2ao + a2) (1) 

and A stands for the regular part of physical scattering amplitude for tt+tt" tt^tt*^ at threshold, calculated at 
the leading order in isospin breaking, and in all orders in chiral expansion. The strong energy-level shift of the 
7r~p atom is also expressed solely in terms of the regular part of physical 'K~p — > 7r~p scattering amplitude at 
threshold, calculated at the leading order in isospin breaking. The relations between the observablcs of hadronic 
atoms and the physical scattering amplitudes obtained by using the nonrelativistic Lagrangian technique, pro- 
vide a desirable generalization of the relations from Ref. [1] that now contains all leading-order isospin-breaking 
corrections. These relations are obtained under a very general assumptions about the mass differences of isotopic 
partners, and do not resort explicitly to the chiral expansion of the corresponding strong amplitudes. In order 
to extract the strong scattering lengths from the physical scattering amplitudes that are actually measured in 
the experiment, one has, however, to use chiral expansion for these amplitudes, and establish such relations in 
a given order in chiral expansion. 

By using chiral expansion, it is demonstrated that the potential scattering theory approach in its present 
form does not account for all isospin-breaking corrections to the strong amplitude. The dominant isospin- 
breaking contributions which are neglected in the potential approach, have been explicitly identified. First of 
all, these are effects of direct quark-photon coupling, that are encoded in the "electromagnetic" counterterms 
of the chiral Lagrangian. Then, there are effects coming from "tuning" of the quark mass so that the common 
mass of the pion triplet in the isospin-symmetric world coincides, by definition, with the charged pion mass. 
We conclude, that for comparison to the results of potential approach, it is necessary to match the latter with 
ChPT in the isospin-broken phase, in order to consistently take into account above effects. 

[1] S. Deser, M. L. Goldberger, K. Baumann, and W. Thirring, Phys. Rev. 96 (1954) 774. 
[2] A. Gall, J. Gasser, V. E. Lyubovitskij, and A. Rusetsky, Phys. Lett. B 462 (1999) 335. 
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with a high accuracy 

L.Afanasyev, M.Jabitski, A.Tarasov, and O.Voskresenskaya 
Laboratory of Nuclear Problems 
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The interaction of 7r"'"7r~-atom (v427r) with matter is of great importance for the DIRAC experiment [1] as 
the y4.27r breakup (ionization) in such interactions is exploited to observe A2tz and to measure its lifetime. In 
the experiment the ratio of the number of 7r+7r~-pair from the ^27r breakup inside a target to the number of 
produced A2Tr (called the probability of breakup) will be measured. The measurement of the ^2Tr lifetime 
is based on the comparison of this experimental value with its calculated dependency on the lifetime. So the 
accuracy of this calculation is essential for the extraction of the lifetime. 

In the DIRAC proposal the probability of ^27r breakup was calculated by evaluating of a set of differential 
equations for populations of states [2]. These equations describe the evolution of the populations, when 
^27r traverses the target, for all atomic states with the principal quantum number less than some limit. The 
initial condition for this set of equations is given by the probability of the A2Tr production in various quantum 
states [3]. As a solution we can get three probabilities for the case behind the target: 1) the summed population 
Pdsc of all discrete states considered in this set; 2) an estimate for the summed population of all other discrete 
states Ptaii; 3) the probability of A2,r annihilation Panh- The remainder is the probability of the A2Tr breakup 
Pbr: 

Phr — 1 Pksc -ftail Panh • 

The accuracy of this procedure is mainly defined by the accuracy of the cross sections for the ^2^ interaction 
and target atoms. Now approaches arc available to calculate these cross sections in the Glauber approximation 
[4] which takes into account all multi-photon exchanges instead of only the single-photon exchange in the first 
Born approximation as used before. This could provide an accuracy in the cross sections at the level less than 
1% and hence almost the same accuracy for Pbr- 

To estimate a limit in the accuracy of the Pbr calculation based on probabilities we have repeated all our 
calculation of ref.[2] using the parabolic basis for the ^27r states description instead of the spherical one as in 
[2]. In both cases we neglect interference effects in the A2Tr description but its contribution should be different 
for another basis. So, the difference in Pbr obtained with two bases shows a principal limit in accuracy for this 
method. In the following table the results of the calculation are shown for the A2Tr momentum 4.7 GeV/c, the 
lifetime 3.0 • 10~^^ s and target thicknesses equivalent in multiple scattering to 30 /im Ta. 





Z 


psph 
-^br 


ppar 
-^br 


APbr/Pbr 


At/t 


Al 


13 


0.223 


0.229 


0.0292 


0.079 


Ti 


22 


0.326 


0.330 


0.0125 


0.030 


Fe 


26 


0.435 


0.438 


0.0069 


0.018 


Ni 


28 


0.470 


0.473 


0.0059 


0.017 


Mo 


42 


0.540 


0.543 


0.0044 


0.015 


Ta 


73 


0.671 


0.673 


0.0026 


0.015 


Pt 


78 


0.704 


0.706 


0.0022 


0.017 



The value of At/t could be considered as the bias of the method. To go lower than 1% we need to consider 
the process of the A2Tr passage through the target in terms of amplitudes instead of probabilities. 

[1] B. Adeva ct al., Lifetime measurement o/tt+tt" atom,s to test low energy QCD predictions, (Proposal to the 

SPSLC, CERN/SPSLC 95 1, SPSLC/P 284, Geneva 1995). 

[2] L.G. Afanasyev and A.V. Tarasov, Yad. Fiz. 59 (1996) 2212. 

[3] L.L. Nemenov, Yad.Fiz. 41 (1985) 980. 

[4] L. Afanasyev, A. Tarasov and O. Voskresenskaya J.Phys.G 25 (1999) B7-10. 
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Pionium interacting with matter: 
I. Formalism for coherent interaction 
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The experiment DIRAC, currently being performed at CERN [1], aims at measuring the lifetime for pionium, 
i.e. 7r+7r~-atoms, in its ground state to better than 10% accuracy thus providing a crucial test for chiral 
perturbation theory. Attaining the required precision hinges on an accurate calculation of all electromagnetic 
processes, competing with the strong interaction, to a precision of a few percent. Therefore we have applied the 
semi-classical formalism (SCA) to calculate the electromagnetic excitation and ionization, i.e. breakup, of the 
pionium to such high precision. 

The SCA has been established [2-4] as a powerful tool to investigate such excitation processes even at 
relativistic energies. In this framework we describe an arbitrary collision process in terms of a trajectory 
associated with the projectile and characterized by an impact parameter b. The amplitude for the transition — 
excitation or ionization — of the target system from an initial state |i) to a final state |/) is given in first order 
of the interaction by the matrix element 



dt{f\Hi^t{Rb{mi), 



(1) 



where Rb{t) describes the projectile's trajectory specified by the impact parameter in the target's rest frame 
(see figure). 



Figure: Semi-classical picture of a projectile with charge Zp 
moving on a trajectory with impact parameter b past a tar- 
get atom A consisting of particles of opposite charge. Kbit) 
describes the projectile's trajectory specified by the impact 
parameter in the target's rest frame. 





X 


Zp 


x' 

V 




t = — oo 


b 




O' 


t = +00 


A^ 









In our case, in the rest-frame of the pionium, the target material's ion can be treated as a classical particle 
moving on a straight- line trajectory Rb{t) = (6, 0,/3cf) at nearly the speed of light, while the pionium at the 
origin is treated quantum mechanically. As we are only interested in pions forming atom-like complexes, their 
relative velocity must be small (of order v^^/c w a). Hence non-relativistic hydrogenic wave fimctions are 
perfectly appropriate for the initial and final states \i) and |/) of the pionium. On the other hand, the complex 
charge distribution of the target atoms is taken into account by including a screening function in the scalar 
potential experienced by the pionium. In our calculations, we use screening functions that reproduce exactly 
expectation values of powers of the radial variable obtained with full Dirac-Hartree-Fock-Slater wave functions 
for the heavy ion or atom [5, 6]. Therefore by neglecting magnetic terms (estimated to contribute no more than 
0.4%), the interaction Hamiltonian Hint reduces to the scalar (Licnard-Wiechert-) potential between the ion 
with charge Zp and the pionium. In the pionium's rest frame we thus have: 



$(f,i) = 



Zpe X - 

27r2 ^ 

k = l 



At 



exp 



(f - Rb{t)) 



(2) 



where the screening parameters Ak and are taken from [6]. Performing the integration over coordinate space 
implied in the matrix element (1) with the two pions positioned at ±r/2 from the pionium's center-of-mass, we 
then obtain the impact parameter dependent transition amplitude in first order of the scalar interaction: 



= -^^4^(2^/ + l)(24 + l) (-1)™/ 
X [l-(-l)^] 



-ruf m rrii 



N 

E- 



s ds 





Be^jb, go, s) 



(3) 
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with go = {Ef —Ei)/{f3Uc), the straight-line trajectory factor B^rn{f>, (ZO) s) [7], and the radial form factors Fj-{k), 
depending on the boTind state or continuum wave functions of the pionium, as appropriate. These form factors 
can easily be evaluated using standard methods as described thoroughly in [8] and [9]. 

Integrating the squared amplitude over the impact parameter b yields the inelastic cross section <Tj^^ for the 
transition between the states |i) and |/) and the total cross section is then obtained by summing over bound 
and integrating over continuum final states — or by using the completeness relation for the set of final states: 

<i = ^ 4V = (^) ' £ s,s[l- FUs)] 

The determination of transition probabilities and cross sections thus reduces to the accurate and fast calcu- 
lation of the radial form factors. 

In [8] we have shown the suitability of our approach under quite general conditions. We have calculated 
the total inelastic cross section in dependence of the kinetic energy of the heavy ion (in the rest frame of the 
pionium) for different target ions. Hereby we have shown that the total cross section, divided by the dominating 
(but trivial) factor 1//?^ in (4), is essentially constant in the energy range of interest to experiment DIRAC, i.e., 
between 2 GeV and 10 GeV. Furthermore we have calculated cross sections for various transitions and break-up 
processes also in the screened Coulomb field of a nuclciis both analytically as well as numerically. 

Our on-going developments now concentrate on including magnetic interaction terms, as well as all higher- 
order perturbation contributions. 



JV 

E 

.fe=i 



(4) 



[1] L.L. Nemenov et al.. Proposal to the SPSLC: Lifetime measurement o/tt+tt -atoms to test low-energy QCD 
predictions, CERN/SPSLC 95-1, SPSLC/P 284. 

[2] J. Bang and J.M. Hansteen, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 31 (1959) 13. 
[3] L. Kocbach, Z. Phys. A 279 (1976) 233. 

[4] D. Trautmann and F. Rosel, Nud. Instr. Meth. 169 (1980) 259. 

[5] G. Mohere, Z. Naturforsch. 2a (1947) 133. 

[6] F. Salvat, J.D. Martinez, R. Mayol and J. Parellada, Phys. Rev. A 36 (1987) 467. 
[7] P.A. Amundsen and K. Ashamar, J. Phys. B 14 (1981) 4047. 

[8] Z. Halabuka, T.A. Heim, K. Henckcn, D. Trautmann, and R.D. VioUier, Nud. Phys. B554 (1999) 86. 
[9] D. Trautmann, G. Baur and F. Rosel, J. Phys. B 16 (1983) 3005. 
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The atomic structure was included in the previous calculation only via the (elastic) form factor F{q'^) [1]. Of 
course the (electrons of the) atoms can be excited as well. A simple estimate is that these incoherent/inelastic 
processes are of the order 1/Z compared to the elastic/coherent one; therefore, in order to achieve the desired 
1% accuracy needed for the experiment, a more precise calculation is needed. 

We derive the cross section for the incoherent contribution starting from PWBA (Plane Wave Born Ap- 
proximation) [2]. The electromagnetic transitions of the atoms and the Pionium arc given through the elec- 
tromagnetic tensor W^"^, which due to gauge invariance and current conservation, only depends on two scalar 
structure functions W'^^\q-^ , qP) and W^^'' {q"^ , qP). This is of course well known, e.g., from electron scattering 
[3]. Taking into account only the scalar interaction and neglecting the magnetic terms, only W^"^^ is important 
and is proportional to the square of the charge transition. One advantage of this approach is that even though 
the relative motion of atom and Pionimri is relativistic, the structure functions W can be calculated in the 
individual rest frames and therefore nonrelativistic approximations can be made. 

We neglect recoil corrections and find the total cross section to be 



-I 



duj dA d q± 



4a^ W2n{^,q^)W2A{A,q^) 



with u) the excitation energy of the Pionium (in its rest frame) and A the one for the atom. The momentum 
square of the photon is given by 



2ojA 



The relative Lorentzfactor is given by 7, /3 = v/c. 

We discuss approximations made to find closed expressions for the total cross section. In them the w and A 
dependence in q^ are replaced by average values. We also find the cross section in the elastic case to be identical 
to the one from the SCA calculation, sec previous contribution. 

To get the total incoherent contribution, we sum over all possible excitations of the atom. This so-called 
inelastic scattering function Sine is then given by 



Sinc{k) = j dAW2{A,k^) = \Fxa{k)\ 



with Fxo the transition form factor. For the calculation we make use of mean-field wave functions calculated 
within a self consistent Dirac-Hartree-Fock-Slater calculation [4] . In terms of the single electron wave functions 
Sine is then given by 

z z 

The often used "no correlation limit" for Sine is found to be not accurate enough. This is due to the fact 
that at the small g^, important in our case, atoms are mainly excited and not ionized; therefore Pauli blocking, 
which is neglected in the no-correlation limit, is important. 

Results of calculations within this formalism are presented in the next contribution. 

[1] Z. Halabuka, T. A. Heim, K. Hencken, D. Trautmann, and R. D. VioUier, Nucl. Phys. B554 (1999) 86. 
[2] T. A. Heim, K. Hencken, D. Trautmann, G. Baur, in preparation (1999). 
[3] T. deForest and J. D. Walecka, Adv. Phys. 15 (1966) 1. 

[4] F. Salvat, J. D. Martinez, R. Mayol, and J. Parellada, Phys. Rev. A 36 (1987) 467. 
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We present results of a detailed investigation of the target-elastic and target-inelastic electromagnetic cross 
sections for pionium scattering off various target materials. Within the closure approximation [1] the total 
electromagnetic cross sections from a given pionium initial state i can be written as 



CTtot = 1677 



)2 fc 
L 



dq± q± 



$^(fc) 



(1) 



where the momentum transfer variables s and k refer to the rest frame of the pionium and the atom, respectively. 
/jP; denotes the monopolc formfactor for the pionium in the state i. For the target-clastic total cross section 
the atomic structure is given by = {Z — Foo(A;))^, where Fqo denotes the coherent formfactor of the 

electronic ground state orbitals only, as the contribution from the (point-like) nucleus has been separated out. 
Wc calculate the atomic formfactor in the framework of the Dirac-Hartrce-Fock-Slatcr model, i.e., using orbitals 
obtained from the numerical solution of the Dirac (or Schrodinger) equation for each occupied orbital. However, 
the coherent formfactor is very well approximated by simple analytical expressions as given by Moliere or Salvat 
et al. [2]. 

In the case of the target- inelastic total cross section, the atomic structure is contained in <I>a(A:) = Si-ac{k) as 
given in the preceding contribution. We determine S-mc using the same orbitals as for the coherent formfactor. In 
a simple "no-corrclation limit" , 5inc(fc) = Z — [_foo(fc)]^/^. This approximation is not sufficiently accurate, as it 
completely neglects Pauli blocking and thus over-estimates the cross section. Our approach with Dirac-Hartree- 
Fock orbitals correctly accounts for Pauli blocking in S'inc(fc)- Moreover it affords determining the contributions 
to the pionium cross sections for each atomic shell individually. This information cannot be extracted from 
tabulated values of formfactors and scattering functions [3] where only the combined contributions of all atomic 
shells are given. 

The diagram on the right shows the integrand of (1), multiplied with the integration variable (thus 
representing more clearly the relative magnitude of the individual contributions to the integral on a logarithmic 
scale for q_i_). 

The target material is Ti {Z = 22), pionium is initially 
in its ground state. Note that the choice of atomic 
excitation energy A affects the integrand only at q±_ 
smaller than the range of dominating contributions. 
The main contributions to the integral come from a 
complex interplay of photon propagator, atomic, and 
pionium structure at 1 < < 100 a.u. Note also that 
the outer shells contribute much more to the target- 
inelastic cross section than the inner shells. Not only 
are there more electrons in outer shells than in inner 
shells, but the contribution of each individual target 
electron is approximately proportional to its princi- 
pal quantum number. Finally, note that the target- 
inelastic cross section is significantly larger than the 
target-elastic divided by Z. Thus neither the "no- 
correlation" limit, nor the simple scaling approxima- 
tion (Tine ~ <^coh/Z arc accurate enough. 

[1] Z. Halabuka, T.A. Heim, K. Hencken, D. Trautmann, and R.D. VioUier, Nucl. Phys. B554 (1999) 86. 
[2] F. Salvat et al., Phys. Rev. A36 (1987) 467; G. Moliere, Z. Naturforsch. 2a (1947) 133. 
[3] J.H. Hubbell et al., J. Phys. Chem. Ref. Data 4 (1975) 471; ibid. 6 (1977) 615; J.H. Hubbell and I. 0verb0, 
ibid. 8 (1979) 69. 
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Another exotic relativistic atom: antihydrogen 
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It is certainly difficult to do experiments with fast neutral atoms. Recently, antihydrogen has been produced 
and detected at LEAR/CERN [1] and Fermilab [2]. Cross-sections for the process 

Z Z ^-Tl^ + e- (1) 

are calculated [3] , see also figure below. The formalism is very similar to the one for calculating the breakup of 
Pionium in matter, see also these miniproceedings. Lamb-shift measurements of hydrogen [4] and antihydrogen 
in flight, which are planned at Fermilab are briefly discussed, see Refs. [4,5,6]. 

Maybe related experimental techniques could also be relevant for the measurement of the Lamb-shift in 
Pionium. Of course, the properties of these atoms have very different scales (lifetime, size, mass, ...). This 
2s-2p energy splitting is an interesting quantity since the scattering lengths ao and 02 are involved in a diflferent 
combination as compared to the one, which enters in the Pionium lifetime [7]. 
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Fig. 1: Cross sections for the reaction eq. 1 for Z = \. The Lorentz factor 7 of the incoming antiproton in 
the Fermilab experiment [2] was in the range of 5.6-6.7, corresponding to an antiproton momentum range of 
5200-6200 MeV/c. A cross section of ais+2s = (1.12 ± 0.14 ± 0.09)pb was found. (The cross section cjns for 
capture into the s-state with principal quantum number n scales as 1/n^.) 

[1] G. Baur at al., Phys. Lett. B368 (1996) 251. 

[2] G. Blanford et al., Phys. Rev. Lett. 80 (1998) 3037. 

[3] H. Meier, Z. Halabuka, K. Hencken, D. Trautmann, and G. Baur, Eur. Phys. J. C5 (1998) 287. 
[4] V. V. Parkhomchuk, Hyp. Int. 44 (1988) 315. 

[5] C. T. Munger, S. J. Brodsky, I. Schmidt, Phys. Rev. D49 (1994) 3228. 

[6] G. Blanford et al., Phys. Rev. D57 (1998) 6649. 

[7] L. L. Nemenov et al., CERN/SPSLC 95-1, SPSLC/P 284, 1999. 
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Multiphoton exchange in pionium-matter interaction 



J. Cugnon^ and D. Vautherin^ 
^Physics Department B5, University of Liege, B-4000 Sart Tilman Liege 1, Belgium 
^LPTPE, Universite Pierre et Marie Curie, 4, Place Jussieu, case 127 F-75252 PARIS CEDEX 05, France 



The interaction of a high velocity pionium atom colliding with an ordinary atom is investigated by studying 
the evolution of the internal wave function driven by the time-dependent Coulomb field generated by the ordinary 
atom. This evolution is governed by the following time-dependent equation: 

where Hq is the intrinsic pionium Hamiltonian and where V{t) is the Coulomb interaction between the pionium 
and the external atom. Using the expansion of ip in terms of the eigenstates tpi of Hq, i.e. 



(2) 



one obtains a set of coupled linear differential equations for the coefficients c^: 

3 

Here, we solve this set of equations in a subset of bound pionium states. We assume that the atom is 
moving along a straight line with a constant velocity, in the rest frame of the pionium. We adopt the dipole 
approximation, i.e. that the Coulomb interaction is reduced to the interaction between the electric dipole of the 
pionium and the electric field generated by the atom. For the latter, the Moliere form factor is used. Equations 
(3) are integrated with proper initial conditions Cj{t = 0) = 5ij, corresponding to the pionium being in state 
tfji initially. The cross section for the excitation from state to state is obtained by integrating over the 
impact parameter h 

/•OO 

Gi^j = dbb \cj{t oo)f. (4) 
Jo 

Keeping the initial values Cj = 6ij in the r.h.s. of the equations (3) yields the usual Born approximation [1]. 
We also show that expression (4) with Born approximation is equivalent to the field theoretical expression for one 
photon exchange, as obtained in ref.[2], when recoil correction is neglected, provided a one-to-one correspondance 
is introduced between the impact parameter and the tranverse momentum transfer. This allowed us also to 
correct the dipole approximation used in V{t) at small impact parameter (details can be found in rcf. [3]). The 
difference between the exact solution of Eq. (3) and the Born approximation can thus be interpreted as due to 
multiphoton exchanges. 

This contribution has been evaluated numerically for a few illustrative cases. The depcndance upon the 
impact parameter, the atomic number of the external atom and the incident energy is investigated. For the 
most intense transition IS — > 2P, the effect of multiphoton exchange is a reduction of the transition cross 
section. This reduction is effective for reduced impact parameter b (= b/a, a being the pionium Bohr radius) 
lying between ~ 0.5 and ~ 2. Let us mention that the bulk of the cross section involves b between ~ 0.5 and 
5 and that the dipole approximation is only effective for b less than unity. 



[1] J.-L. Basdevant, "Mecanique Quantique", Ellipses, 1986 
[2] S. Mrowczyhski, Phys. Rev. D36 (1987) 1529 

[3] J. Cugnon and D. Vautherin, Proc. of "Hadronic Atoms and Positronium in the Standard Model", ed. by 
M.A. Ivanov et al, 1999, p. 128 
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The DEAR experiment at DA$NE 

Carlo Guaraldo 
(on behalf of the DEAR collaboration) 
INFN-Laboratori Nazionali di Frascati, CP. 13, Via E.Fermi 40, 1-00044 Frascati, Italy 

A new era in the field of low energy kaon physics will begin with the start of the new DA$NE low-energy, 
high-luminosity e+e" collider. The objective of the DEAR ( DA$NE Exotic Atom Research) experiment is 
the determination of the isospin dependent KN scattering lengths via the measurement of the strong interaction 
shifts and widths of the kaonic hydrogen and kaonic deuterium i^T-series lines [1]. In practice, in the case of 
kaonic hydrogen, the Kq, line at ^ 6.15 keV is the most important one for the determination of the shift and 
width of the Is level. 

The challenging goal DEAR is aiming at is to measure the line shift in hydrogen to a precision of 1 %, 
the width to a precision of a few percent, and kaonic deuterium for the first time. 

The main challenge of DEAR is to detect a very weak x-ray signal out of the background of an c+c~ collider 
with high statistics and high precision. Only recently these X rays have been observed for the first time (KpX 
at KEK [2]) with an overall statistics of only 114 counts in the peak. 

To overcome the background problem Charge-Coupled Devices (CCDs) arc used as x-ray detectors. CCDs 
are characterized by a good energy resolution and by an unprecedented background rejection capability [3] . 

A kaonic atom is formed when a negative kaon enters a target, loses its kinetic energy through ionization 
and excitation of the molecules of the medium and eventually is captured, replacing the electron, in an excited 
orbit. Various cascade processes deexcite this kaonic atom to the ground state. 

When a kaon reaches low-n states with small angular momentum, it is absorbed through a strong interaction 
with the nucleus. This strong interaction causes a shift in the energies of the low- lying levels from their purely 
electromagnetic values, whilst the finite lifetime of the state turns out in an increase in the observed level width. 

The shift e and the width T of the Is state of kaonic hydrogen are related in a fairly model-independent 
way to the real and imaginary part of the complex s-wave scattering length, a^-pi e + = 20? ii^a^-p = 
(412 eVfm~^) • a^-p (Deser-Trueman formula [4]), where a is the fine structure constant and /x the reduced mass 
of the K~p system. A similar relation applies to the case of kaonic deuterium and the corresponding scattering 
length, a^-d- 

These observable scattering lengths are related to the isospin dependent scattering lengths ag and ai in the 
following way: ax-p — ^(ao + ai) , o-K-d = \ {^ rnN+mK/2 ) ("0 ^^i) + C .\n the case of deuterium, the 
first term represents the lowest-order impulse approximation, in which the kaon scatters on each "free" nucleon. 
The second term, C, contains all high-order contributions, including three-body effects. This term can be larger 
than the first term, thus the extraction of the two isospin scattering lengths from an observation of the kaonic 
hydrogen and kaonic deuterium scattering lengths requires a dedicated analysis. An accurate determination 
of the N isospin dependent scattering kuigtlis will place strong constraints on low energy K~ N dynamics, 
which in turn constraints the SU(3) description of chiral symmetry breaking. 

Crucial information about the nature of chiral symmetry breaking, and to what extent chiral symmetry is 
broken, is provided by the calculation of the meson-nucleon sigma terms. The meson-nucleon sigma terms are 
defined as the expectation value of a double commutator of the chiral symmetry breaking part of the strong- 
interaction Hamiltonian. The low energy theorem relates the sigma terms to the meson-nucleon scattering 
amplitude. A phenomenological procedure starting from the experimental amplitudes is then used to determine 
the sigma terms and therefore measure chiral symmetry breaking. According to an evaluation based on the 
uncertainties in the phenomenological procedure the sigma terms can be extracted at the level of 20 %, by 
combining the precision measurement to be performed by DEAR with the bulk of most recent analyses of 
low energy K^N scattering data. The sigma terms are also important inputs for the determination of the 
strangeness content of the proton. 

Presently, the DEAR experiment is installed on DA$NE and ready to start data taking. 

[1] S. Bianco et al., LNF-Preprint, LNF-98/039 (P) (1998); Rivista del Nuovo Cimento, November 1999. 
[2] T.M. Ito et al., Phys. Rev. A 58 (1998) 2366. 
[3] J.-P. Egger et al., Part. World 3 (1993) 139. 

[4] S. Deser et al, Phys. Rev. 96 (1954) 774; T.L. Truemann, Nucl. Phys. 26 (1961) 57. 
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Pionic hydrogen: status and outlook. 



Leopold M. Simons (for experiment R-98-01.1 at PSI) 
Paul Scherrer Institute, 5232 Villigen PSI, Switzerland 

The measurement of the strong interaction shift and width of the ground state of the pionic hydrogen atom 
allows to determine two different linear combinations of the two isospin separated s-wave scattering lengths of 
the pion nucleon system. Past experiments at PSI measured the 3-1 transition energy in pionic hydrogen and 
deuterium [1,2,3]. The measurement combined a highly efficient stopping arrangement (cyclotron trap) and a 
sophisticated cylindrically bent Bragg crystal set-up with newly developed CCD detectors. The shift was finally 
determined with an accuracy of better than 10~^. The error in the determination of the width of the ground 
state level in hydrogen is a factor of almost one order of magnitude worse than the error in the shift value, 
which inhibits a determination of both isospin separated scattering lengths on the percent level from the pionic 
hydrogen measurement alone. The situation can be improved in a combined analysis with the pionic deuterium 
shift measurement, but then one has to rely on a theoretical understanding of the pionic deuterium system. 

A new proposal at PSI aims at a direct determination of the width of the ground state of pionic hydrogen 
on the level of 10~^ [4]. In order to achieve this the experiment should meet several requirements as are 
a higher number of detected X-rays and a thorough understanding of the response function of the crystal 
spectrometer. Most importantly it should also be able to disentagle the strong interaction broadening from 
other line broadening effects as the Doppler effect. Pionic hydrogen atoms may gain kinetic energy by converting 
binding energy difference into kinetic energy, a process known as Coiilomb dcexcitation. A recent experimental 
investigation of this effect took place at PSI. The distortion of neutron time of flight spectra from the charge 
exchange reaction has been measured in liquid hydrogen and in gas at a pressure of 40 bar (room temperature) 
[5] . The measured spectra clearly show the effect of different transition steps converted into kinetic energy. 

The statistics problem was solved by stopping more pious in a newly built cyclotron trap and by increasing 
the efSciency of the crystal spectrometer by using spherically bent crystals as well as bigger CCD detectors. 
In this way a factor of more than one order of magnitude was gained in count rate. In addition the peak to 
background ration as well as the resolution of the spectrometer had been improved. First experiments with 
pionic deuterium show the high quality of the data achieved by these measures [6]. 

For a better determination of the response function an ECR (Electron Cyclotron Resonance) source is 
presently being set up which produces X-rays from hydrogen-like ions with high intensity. Their line widths 
are expected to be an order of magnitude narrower than the expected resolution of the spectrometer. Under 
these circumstances a detailed measurement of the response function and even an optimization of the crystals 
will be feasible. 

The third problem will be attacked in two steps. First the 2-1 and the 3-1 transitions in pionic hydrogen 
will be measured at about three different pressures. The influence of the Doppler effect on the line shape will be 
different for the different transitions and pressures whereas the strong interaction width will remain the same. 
Monte Carlo simulations have shown that a suitable fitting procedure will allow to reach an accuracy in the 
width determination of 3 x 10~^. 

In a second step a simultaneous mesurement of muonic and pionic hydrogen X-rays will be performed. The 
muonic hydrogen X-rays are affected by Doppler effect in a similar way as their pionic counterparts. In contrast 
to them they are not broadened by strong interaction, however, and therefore offer a good chance to determine 
the influence of the Doppler effect [7] . A further check of the theory of the cascade processes will be carried out 
by measuring the velocity of the system at the instant of the charge exchange process by the neutron time of 
flight method. The pionic X-ray measurements will start end of year 2000 and first results can be expected in 
the year after. 

[1] D. Sigg et al., Nucl. Phys. A 609 (1996) 269. 

[2] D. Chatellard et al., Nucl. Phys. A625 (1997) 855. 

[3] H.-Ch. Schroder et al., to be published in Phys. Lett. B 

[4] Proposal R-98-01.1 at PSI. 

[5] A. Badertscher et al., Phys. Lett. B392 (1997) 278. 
[6] P. Hauser et al. Phys. Rev. C58 (1998) R1869. 
[7] V. E. Markushin, private communication. 
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Pion— kaon scattering 



Ulf-G. Meifiner 
FZ Jiilich, IKP(Th), D-52425 Jiilich, Germany 



Pion-kaon scattering has been considered in refs.[l,2] in the framework of three flavor chiral perturbation theory 
to one loop accuracy. To that order, all low energy constants have been determined from other reactions (and 
invoking large Nc arguments), so that phase shifts, threshold and subthreshold parameters can be predicted. 
Despite the rather large threshold energy, ^/s^^^^. = Mk + « 633 MeV, the one loop corrections to the S-wave 
scattering lengths are about 25%. The uncertainties due to the various parameters at next-to-leading order 
are also discussed in [1] . Other approaches and extensions (explicit resonances [3] , heavy kaon approach [4] , the 
inverse amplitude method [5] extended to coupled channels [6]) all lead to very similar threshold parameters. 
The experimental values scatter over large ranges, see e.g. the figure in [2]. Also, two groups have performed 
a dispersion theoretical analysis [7,8]. While a^^"^ and a\^'^ agree with the chiral predictions, the dispersion 

1 /2 

relation result for is sizeably larger. A new dispersive analysis including also the SLAG data from 1988 is 
in progress [9]. In that context, we have also derived a sum rule. 



Even if one assumes that Im F~ is entirely given by the if* (892) and approximating the latter by a (5-function, 
one obtains FkF^^ « (127 MeV)^ not far off the empirical value of (102 Me V)^. Other work in progress contains 
a one loop analysis of 'k~ tt^K^^ [10], which is needed to calculate the width of nK atoms. Finally, an 
analysis of ttK scattering in GCHPT is also available [11]. 
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decays at DA$NE 



Patrizia de Simone 

for the KLOE Collaboration 
Laboratori Nazionali di Frascati dell'INFN, Frascati, Italy 

The capabilities of the KLOE detector [1] , and the over constrained kinematic of the (f>'s decays at DA$NE 
[2] will give the opportunity to improve our knowledge on the K14 decays. 

The experimental results obtained so far are dominated by the work of L. Rosselet et al. [3] which measures 
the TT+TT" phase shifts with ~ 30000 K^^ events collected. The phase shift difference {5° — Si) was determined 
in five bins of the invariant mass of the dipion system, and the related isoscalar S-wave scattering length 
was a° = 0.26 ± 0.05. This result must be compared with the CHPT prediction a° = 0.20 ± 0.01 [4]. The 
apparent discrepancy between the CHPT prediction and the L. Rosselet experiment may be interpreted as a 
manifestation of an Timisual low value of the quark-antiquark condensate: a new measurement of the isoscalar 
S-wave scattering length a° with a substantially smaller error can determine if this discrepancy is statistically 
significant. 

Two methods have been proposed to measure the isoscalar S-wave scattering length a°: the Pais Trei- 
man method [5], and the Maximum Likelihood Method (MLM) [6]. The MLM requires to do approxima- 
tions/assumptions on the form factors, while the Pais Trciman method is model independent. The two methods 
are complementary and both will be considered at KLOE. Moreover, since the Pais Treiman method was devised 
for a detector with a uniform efficiency over the whole phase space such as we expect KLOE to be, this will 
be the first time that this method will be applied (a variant of the Pais Treiman method has been used by L. 
Rosselet et al. [3]). 

The statistical errors on the Kf^ parameters has been evaluated using the MLM method [6], in particular, 
concerning the tttt phase shifts, the error on the isoscalar S-wave scattering length has been estimated using 
the J.L. Basdevant et al. parametrisation [7]. For 300000 K^^ events the estimated error is 6a° = 0.01. This 
estimate is purely statistical and apply to a "perfect" detector, i.e., one which covers the whole phase space 
with unity efficiency everywhere. This is close to being true for KLOE, that it is a hermetic detector operating 
at DA$NE producing self-tagging pairs, with high reconstruction efficiency of neutral and charged low 
energy particles that will be controlled at the level of ~ 10~^ 10~^. 

In one year of running at DA<I>NE at the target luminosity of £0 = 5 x 10"^^ cm~^s"^, given that the cross 
section for e+e" ^ ^ at the (p resonance peak is ~ 3.8^b, the number of (j) decays into charged kaons will be: 

iV(0^is:+i4:-)~9.4x 10^ 

and the number of K^^ decays are: 

N{K^^) = 1.88 X 10^° X 3.9 x 10"'^ x ^ x ^ ~ 4.2 x 10^ 

including the efficiency e/f± to select the (j) decays in charged kaon events, and the cuts {cminv) the Dalitz 
invariant mass plot. To collect 300000 K^^ decays in one year, the reconstruction efficiency of the K^^ events 
has to be ~ 71%. 

The KLOE detector at the Frascati ^-factory DA$NE is fully operational and has recently started to collect 
its first data. During a short test run period 220 nb^^ of integrated luminosity has been accumulated and has 
been used to study the detector performance and to fully test the reconstruction program. The detector is 
ready to collect and reconstruct data at the expected luminosity. 

[1] J. Lee-Pranzini, The second DA$NE Physics Handbook, ed. L. Maiani et al. (Frascati, 1995) p.761. 

[2] G. Vignola, Proc. Workshop on Physics and Detectors for DA$NE, ed. G. Pancheri et al. (Frascati, 1991) 

p.U. 

[3] L. Rosselet et al., Phys. Rev. D15 (1977) p.574. 

[4] J. Gasser and H. Leutwyler, Phys. Lett. B125 (1983) p.321. 

[5] A. Pais and S.B. Treiman, Phys. Rev. 168 (1968) p.l858. 

[6] M. Baillargeon and P.J. Franzini, The second DA$NE Physics Handbook, ed. L. Maiani et al. (Frascati, 
1995) p.413. 

[7] J.L. Basdevant et al, Nucl. Phys. B72 (1974) p.413. 
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A new measurement of the tt+tt e^u decay. 

Stefan Pislak for the E865 collaboration 
Yale University 



The measurement was performed by experiment E865 at the Brookhaven Alternating Gradient Synchrotron 
employing an apparatus that has been described in previous publications [1]. The selected A'e4 candidates 
were required to have an unambiguous identification of the to suppress background originating from 
decays with a 7r° (mostly Kt^2i ^eSj and giving an e+e^ pair. In order to reject the — > Tr^'^Tr+Tr^ 

{Kt) background, we require an unambiguous identification of the e+ and a kaon, reconstructed from the three 
charged decay products, which does not track back to the target. The basis of the rejection of accidental 
background is the requirement to have three trajectories consistent with having come from a common vertex, 
and a timing spread between the trajectories consistent with the resolution. Our final signal sample contains 
437,000 K^j^ candidates including 2% background events. These data represent a more than tenfold increase 
in statistics compared with previous experiments [2]. Fig. ^) shows the reconstructed tttt mass in comparison 
with the Monte Carlo simulation, one of many control plots, demonstrating our good understanding of the data. 

These data were then used to fit phase shift differences and form factors as a function of the tttt mass. In 
a first analysis, a formalism as it was described by Rosselet [2] was employed. Radiative corrections were not 
included but their calculation is in progress [3] . A preliminary fit of the phase shift differences as a function of the 
TTTT mass is shown in Fig. ^p). An estimation of the phase shifts and form factors based on the parameterization 
proposed by Amoros and Bijnens [4] and the determination of the branching ratio is in progress. 
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Figure 3: a) Distribution of the reconstructed tttt mass. The dots represent the data, and the histogram shows 
a Monte Carlo simulation, b) Dependence of the phase shift difference 5 = 5^ — 5\ on the tttt mass. The fits are 
based on the work by Basdevant et al. [5] . 



hep-ex/9907045, to be published in Phys. Rev. Let. 83, 29 November 1999. 
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K£4: decays: a theoretical discussion 



G.Amoros 
Lund University 



To measure the s-wave of the tttt scattering, the experiments need to focus on processes where the inelastic 
channels are supressed. The decay of the pionium system and the _ft'£4 decays are competitive and independent 
ways. This talk will be related with K^^ decays where the supposed precision for the experiment will force 
the calculation to higher orders in CHPT [1]. An overview of the status of the theoretical calculation and the 
influence on the experimental parametrization [2] is given. A first part of the talk is related with a general and 
model independent parametrization for the K(4 decays. Previous study of the prediction for the relevant form 
factors in several models and in the framework of CHPT allows us to consider the linear dependence with the 
kinematical variables as the natural one. This parametrization could decrease the number of variables to fit. In 
the second part of the talk some results about the two-loops calculation is given. Decay constants and masses 
for the complete SU{3) octect are obtained. 

[1] G. Amoros, J. Bijnens and P. Talavera, work in preparation. 
[2] G. Amoros and J. Bijnens, J. Phys. G. 25 (1999) 1607. 
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Chiral phase transitions at zero temperature 



Jan Stern 

I.P.N. , Groupe de Physique Theorique, Universite de Paris-Sud, F91406 Orsay Cedex , France 



It is argued that in QCD , quark loops progressively suppress order parameters of chiral symmetry breaking 
as the number Nf of light flavours increases. This could result into a non-trivial chiral phase structure of the 
theory considered as a function of Nf and N^. For Nf close to the first critical point ncrit{Nc), one may expect 
[1] a reduction of the size of the chiral condensate < qq > and an enhancement of its fluctuations leading to 
an important violation of the Zweig rule (i.e. to a breakdown of largeiVc predictions) in the scalar iso-scalar 
channel. Starting from observed correlation between uu and ss, it is possible to estimate the fall-off of the 
condensate < uu >isij as the number of light spectators varies from Nf = 2 to Nf — 2> [2] . Existing data are 
consistent with a strong Nf - dependence : The three- flavour condensate could be suppressed by as much as a 
factor 2 relative to the two-flavor condensate, possibly indicating that the real world {Nf = 2 — 3, TVc = 3) could 
already feel the influence of a nearby phase transition [2] . It is important to make a careful distinction between 
< uu >2 and < uu >3 in all phenomenological analysis. The two-flavour condensate (in appropriate units) is 
measurable in precise low-energy vr — tt scattering experiments and is strongly correlated with the quark mass 
ratio r — 2mg/ (m„ -|- m^) [1] . The main goal of these experiments remains the determination of the two- flavour 
condensate, in order to examine to which extent the expected decrease of the condensate < uu >Nf towards 
the critical point already starts to be visible at Nf = 2. Whatever the experimental answer to this question 
will be, additional experimental information will be necessary to pin down the three-flavour condensate and to 
settle the theoretical issue of a nearby phase transition as a common explanation of the observed large Zweig 
rule violation in the vacuum-channel and of the Nf dependence of the chiral symmetry breaking. This might 
require new experiments involving Kaons ,in which the strange quark would play a more direct role than that 
of a "sea-side-spectator" . 

[1] S. Dcscotes, L. Girlanda, J. Stern: "Paramagnetic Effect of Light Quark Loops on Chiral Symmetry Break- 



ing" , |hep-ph/9910537[ 

[2] B. Moussallam: " Nf Dependence of the Quark Condensate from a Chiral Sum Rule", ^ep-ph/9909292| 
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The two-flavor chiral condensate from low-energy tttt phase-shifts 

Luca Girlanda 

I.P.N. , Groupe de Physique Thcoriquc, Univcrsitc dc Paris-Sud XI, F-91406 Orsay Cedex, France 



Low-energy tttt scattering has been identified as the most promising physical observable for measuring the 
quark condensate [1]. What can actually be measured in tttt experiments is the two- flavor quark condensate, 
i.e. (qq) in the SU(2)xSU(2) chiral limit [2]. Due to the properties of analyticity, unitarity and crossing 
symmetry and to the fact that Goldstone bosons interact weakly at low energy, the amplitude A{s\,t,u) can 
be written, up to and including the two-loop level, in terms of six parameters (a, /?, Ai, A2, A3, A4). This low- 
energy representation holds independently of the size of (qq) . The explicit expression can be found in Ref . [3] . 
a and f3 represents essentially the amplitude and its slope at the symmetrical point s = t = u = 4/3M^. 
They can be directly related to the two-flavor quark condensate measured by the Gell-Mann-Oakes-Renner 
ratio X2 = —2rh\vaim^o{qq)/{F^M'^). In addition to constraints from chiral symmetry, the amplitude satisfies 
also a set of crossing-symmetrical dispersion relations, known as Roy Equations. Using as input experimental 
phase-shifts in the medium- and high-energy region, they allow to determine the low-energy amplitude in terms 
of two subtraction constants, usually identified with the two scalar scattering lengths ag and a§. In fact the two 
subtraction constants can be related to each other, imposing the matching of the amplitude to the experimental 
P-wave, which is well known down to rather low energies, thanks to the p resonance. This results in the so-called 
Morgan-Shaw universal curve, that we quote in the form given by Petersen in 1979, 
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By comparing the two low-energy representations one is able, in principle, to determine 4 out of the 6 parameters 
of the chiral amplitude. A similar program has been undertaken in Ref. [3], in order to fix the parameters 
Ai, . . . , A4. It turns out that the outcome of this analysis is only barely dependent on a and (3, and hence on 
the quark condensate. The slight a- and /3-dependence can be conveniently approximated by linear functions 
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while the experimental uncertainties on the Ai are independent of a and /3, S}^^"^ — 2.2 • 10^^^ 6X^2^ = 0.5 • 10-^, 
JAg'^P = 0.6 • 10~*, SX^^"^ = 0.12 • 10-**. With the parameters Ai, . . . , A4 determined from the Roy dispersion 
relations, the low-energy tttt observables are effectively parametrized only by a and /3, which can be fitted 
to the experimental phase-shifts. However, in order to fully exploit the analytical properties of the amplitude 
contained in the two-loop six-parametrical formula of Ref. [3] , we can treat all the six parameters as fit variables, 
and add the X°^^{a,l3) as additional experimental points to the . Moreover one can insert the constraint of 
the Morgan-Shaw universal curve, which, using the two- loop expressions for the scattering lengths [3], can also 



be expressed in terms of the fit variables (a, /3, Ai 



, A4). Thus we define the x as 



X = 



E 



(5„(a,/3, Ai,A2,A3,A4) -5°^^ 



4 r 



E 



5XT^ 



+ XMorgan-Shaw(Q^i Xi, A2, A3, A4), 



where J^'^p are the experimental phase-shifts, with experimental error cr'^^ . It is straightforward to include 

2§, once the pionium lifetime will be measured. Applying this fit procedure 



to the the combination a[j 
to the old Kf.i Rosselet data [J^'^p = 5Q{En) — 5\{En)\, we find for the 5- wave isoscalar scattering length, 
Oq = 0.268 ± 0.043, and for the combination relevant to the pionium lifetime, Oq — ap = 0.294 ± 0.033. Notice 
that the error on the difference of scattering lengths is smaller than the individual error on Oq, reflecting the 
importance of properly taking into account all the correlations among the fit variables. Using the relationship 
between X2 and {a, (3), known at one-loop level, we get X2 = 0.71 ± 0.13cxp ± 0.05thoo, where the experimental 
and theoretical errors can be considered as mutually uncorrelated. 
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Semi-lcptonic decays are an important source of information on liadronic matrix elements of the Standard 
Model vector and axial currents and, in the case of the e^v^ decay, on low-energy tt — tt scattering. 

In order to meet the accuracy of two-loop calculations and the expected improvements on the experimental side 
in the near future (BNL-E865, KLOE, DIRAC,...), electromagnetic corrections to these processes, or isospin 
breaking contributions in general, have to be estimated in a reliable way. 

The first step towards this aim has been proposed in Ref . [1] . In order to address the question of electromag- 
netic corrections in semi-leptonic processes on a systematical basis, the particle spectrum of chiral perturbation 
theory with virtual photons has been enlarged by including also the light leptons. The presence of radiative 
corrections affects the semi-leptonic amplitudes in two ways. First, the factorization property between the 
hadronic and the leptonic weak currents is lost, due to loop graphs involving both virtual photons and virtual 
leptons. Second, these additional loop graphs generate in general new divergences, which are taken care of 
neither by the Gasser-Leutwyler low-energy constants Li [2], nor by Urech's counterterms Ki [3]. The structure 
of these divergences have been analysed at one loop, using super-heat-kernel methods [4], and a complete list 
of the corresponding counterterms has been given in [1]. As an illustration of this framework, an explicit cal- 
culation of the decay rates for tt — > £v{'y) and K ^i^i'j) has also been provided. The extensions to the other 
semi-leptonic processes, and to the pion (3— decay, are in progress. 

A numerical estimate of the corresponding decay rates eventually requires a determination of the countert- 
erms involved in each amplitude. This may be achieved following methods originally suggested in Ref. [5], 
combined with large- iVp techniques [6]. From this point of view, the recent study [7] of the decay of pseu- 
doscalar mesons into pairs of charged leptons, which, as far as some aspects of this question are concerned, may 
be considered as a "weak SU(2) rotation" of the tt ^ decay, offers interesting perspectives. 
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Numerical solutions of Roy equations 

Gilberto Colangelo 
Institut fiir Theoretische Physik der Universitat Zurich 
Winterthurerstr. 190, 8057 Zurich 

We solve numerically the Roy equations [1] for the S and P waves of the nn scattering amplitude in the 
low-energy region (below 800 MeV) [2] . The imaginary parts above 800 MeV or of higher waves are taken either 
from experiments (wlica'c; available) or from theoretical modelling. We arc able to solve the equations for any 
S wave 7 = scattering length a^. At fixed Uq solutions can be found only for Uq (the S wave 1 = 2 scattering 
length) within a certain range, called the Universal Band. 

We then confront the solutions wc find for various values of ajj and ag to the available experimental data. 
The corresponding analysis yields ellipses of given confidence level inside the Universal Band. We confirm 
old analyses [3] of the Ke4 data [4] leading to an oq between 0.20 and 0.30 [5]. This machinery is ready to be 
used with forthcoming Ke4 data which should reduce drastically the error on a^. 
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We discuss the use of Non Relativistic Lagrangians for the study of bound state characteristics and in 
particular of hadronic atoms. This study can give important information about strong interactions at low 
energies. Hence it has been the main aim of many experiments performed in the last years and of others, like 
DIRAC at CERN and DEAR at DA$NE, which are running now or planned to start soon. The nature of the 
problem (and mainly the fact that the bound state components momenta are much smaller than their masses) 
suggests the use of non relativistic lagrangians. This approach, adopted in a number of recent publications 
[1-2], has proved to be probably the most efficient way to study these systems in quantum field theories. In this 
connection, we feel that it is useful to discuss general properties of the non relativistic framework in a detailed 
and systematic manner. We refer in particular to the correct matching and regularization of the theory and to 
the problem of ensuring the UV finiteness of the bound state characteristics and their independence on the the 
off-shell behaviour of the Green functions. 

We found it worthwhile to study these problems by using an appropriate model. Our results [3], that will 
appear soon, have proved the possibility of developing a consistent quantum field theoretical framework for 
studying bound states with non relativistic lagrangians. Our model has a massive scalar bound in an external 
coulombic field and interacting with another massless scalar field. It has the advantage of being exactly solvable 
and it presents all the essential features of the physical systems one wants to study, without being affected by 
some inessential complications. We have studied this model at different levels of complexity. 

Neglecting, at first stage, the interaction with the massless scalar, one computes the relativistic wave func- 
tions and energy levels of the heavy scalar in the external field, by looking for the poles of the Fourier transform 
of the two point function. In a non relativistic framework one recovers an effective lagrangian starting from 
the relativistic one and requiring that the S matrix elements are the same in both theories. Then, one can use 
Rayleigh-Schrodinger perturbation theory to evaluate the energy levels. Our lagrangian, written as an expansion 
in inverse powers of the heavy mass, is built in such a way that the non relativistic and relativistic Green func- 
tions are the same and this guarantees that the correct energy levels are reproduced in the effective theory. The 
lagrangians obtained with different techniques give the same S matrix elements on mass shell. We have shown 
that they differ only by terms that can be eliminated by using field redefinitions and that don't contribute to 
the bound state energies. We have also verified that, even in the effective theories, the divergences appearing in 
the bound state calculations can be cancelled, by summing all the contributions of a given order in the coupling 
constant and using the same counterterms needed to make the Green functions finite at the same order. 

Finally, we have added an interaction term between the heavy scalar field and a dynamical massless scalar 
and computed at first order in the coupling constant the energy shift (due to radiative corrections) for the ground 
state of the heavy scalar. To perform the non relativistic calculation one has to face the problem of the breaking 
of the naive power counting rules. In, fact computing the self energy ^ in dimensional regularization one sees 
that the loop contributions fail to reproduce the expansion in powers of the momenta on the basis of which the 
effective lagrangian is built. So one must impose some additional rules, like the so called "multipole expansion" 
that we have adopted, recovering an efficient power counting. With this prescription we have verified that the 
non relativistic energy shift at the first order in the coupling constants reproduces exactly the relativistic one. 

To summarize, it is possible to develop a completely consistent and well defined formalism, based on the use 
of non relativistic lagrangians, to study the hadronic atoms. These techniques already applied successfully to the 
case of pionium [1-2] could be extended also to the analysis of other exotic atoms (like pn^ , pK^ , dir^ ,dK^). 
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^ The problem of self energy calculation and of a correct definition of the mass in the effective theory will be treated in detail in 
a separate work [4] 
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